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We propose a method for implementation of an universal set of one- and two-quantum-bit gates 
for quantum computation in the system of two coupled electrons with constant non-diagonal ex¬ 
change interaction. Suppression of the exchange interaction is offered to implement by all-the-time 
repetition of single spin rotations. Small g-factor difference of electrons allows to address qubits 
and to avoid strong magnetic field pulses. It is shown by means of numerical experiments that 
for implementation of one- and two-qubit operations it is sufficient to change the amplitude of the 
magnetic field within a few Gauss, introducing in a resonance one and then the other electron. To 
find the evolution of the two-qubit system, we use the algorithms of the optimal control theory. 

PACS numbers: 73.21.La, 85.75.-d, 03.67.Lx 


INTRODUCTION 

There have been numerous proposals for implementa¬ 
tion of quantum computation schemes in different mate¬ 
rial realizations, such as photon qubits, trapped atoms 
and ions, nuclear spins in molecules in liquid solutions, 
spin or charge states in quantum dots (QDs) or dopants 
in solids, and superconducting circuits |l|. Each such sys¬ 
tem must satisfy to 5 DiVinchenzo criteria: scalability, 
initialization ability, long coherence lifetime, universal set 
of quantum gates realization, readout ability 0 . Inspired 
by the Loss and DiVincenzo proposal Q, there is a con¬ 
tinuing experimental effort to realize electron spin quan¬ 
tum bits or qubits in semiconductor quantum dots.Q 
There, a single qubit is defined by the two spin states 
of an electron localized in a quantum dot (QD). Cou¬ 
pling between qubits is provided by exchange interac¬ 
tion between electrons in neighboring quantum dots. For 
singie-qubit rotations the electron spin resonance (ESR) 
technique was proposed Q. Two-qubit operations can 
be performed as combinations of the V SWAP opera¬ 
tion (based on the exchange interaction) with single-spin 
rotations. Read-out procedure can be realized via the 
spin-to-charge conversion Q. In light of the foregoing, 
semiconductor QD systems are well suited for quantum 
computing, provided that the decoherence time in the 
system is long enough and there is a possibility of the 
system initialization. Realization of this proposal needs 
a system of electrodes. One type of electrodes provides 
the control of the exchange coupling during two-qubit 
operations, and the second type of electrodes allows to 
address an individual spin for one-qubit operations. In 
subsequent articles, the g-factor engineering or the local 
magnetic field gradient were proposed additionally for 
one-qubit operations 0, [1| . 
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Two-qubit manipulation using electrical gates is actu¬ 
ally quite a challenge. The complexity of this approach 
can be recognized by comparison of two following dates, 
the Loss and DiVincenzo proposal @ has been published 
in 1998, while an experimental implementation of two- 
qubit operation in a double quantum dot has been made 
only in 2011 (see Ref. H). All used electrodes can be 
sources of fluctuations and can lead to undesired errors. 
The solution of this problem is to exclude the gates at 
least partially, for example, to refuse electrodes control¬ 
ling the exchange interaction and to find a way of the 
coupling control without electrodes. 

There are some works where authors proposed quan¬ 
tum computation schemes based on the constant ex¬ 
change coupling. The exchange coupling was eliminated 
by using encoded bits ilH, where two physical qubits 
in the states f J, are unified in one logical qubit 
t/. This configuration provides compensation of inter¬ 
actions with the environment. But two-qubit operations 
demand leaving the interaction-free space. By a strong 
local magnetic field the state fi transforms to ft) and the 
interaction between logical qubits becomes switched on 
Q. After performing the two-qubit operation the qubits 
are driven back to the interaction-free space. One more 
approach for elimination of the exchange interaction is 
based on refocusing pulses that are used in nuclear mag- 
netic resonance (NMR) technique [Sec. 7.7.3 in Ref. [ll|. 
It is applicable if the inter-qubit interaction Hamiltonian 
is diagonal in the computational basis. The repetition 
of refocusing pulses allows to suppress the interaction 
between qubits. Refocusing requires very fast repeated 
switchings with a period much shorter than the elemen¬ 
tary operation time, that is quite difficult to implement 
in practice. Another complexity is locality of refocus¬ 
ing pulses, it requires injecting pulses onto every qubit of 
the quantum system. The use of strong local magnetic 
fields or the refocusing technique requires a more com¬ 
plicated architecture of the quantum computer, and we 
have a goal to find another way to eliminate the coupling 
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between qubits. 

An idea proposed in the present work is close to the 
one of the work of Ozhigov and Fedichkin, [T^ where 
elimination of the exchange interaction by all-the-time 
repetition of one-qubit operations was proposed. We of¬ 
fer repeating single-spin rotations for suppression of the 
exchange interaction, when the latter must be switched 
off. Our approach allows not only elimination of the ex¬ 
change interaction, but also implementation of all basic 
logical operations. 

This paper is organized as follows. In Sec. I we describe 
the intuitive scheme of our proposal, that is represented 
as a mathematical model in Sec. II. In Sec. Ill the ways of 
implementation of basic logical operations are discussed 
in the frame of this model. In Sec. IV the possibil¬ 
ity of one- and two-qubit operations implementation is 
demonstrated by means of numerical experiments. The 
evolution of qubit states is found using the algorithms of 
the optimal control theory [ 1 . Optimization was per¬ 
formed using a combination of nonlocal search methods 
and a local descent by conjugate gradient methods. 


I. MAIN IDEA 

The intuitive scheme of our proposal is the follow¬ 
ing. We take two electrons with a constant exchange 
coupling between them and provide all-the-time rotation 
of their spins by application of the resonance microwave 
frequency radiation. These spin rotations eliminate the 
exchange interaction between electrons. For two-qubit 
gates implementation we switch on the exchange coupling 
taking electrons out of resonance. This can be done by 
simple addition of a small magnetic field B ± 6B{t). The 
most important two-qubit operation CNOT can be done 
through the combination of single spin rotations and the 
VSWAP operation. For the single-spin rotatio ns elec- 
trons are driven to the resonance, and for the \/SWAP 
operation electrons are taken out of the resonance. Sin¬ 
gle spin rotations (when one electron spin of the coupled 
pair rotates, and the second one does not rotate) need the 
difference in electron g-factors. If this difference would 
be not so large (for example, 5g ^ 10“^) one can perform 
quantum operations by a small addition to the magnetic 
field 5B{t) ^ 1 G, that can be produced by conventional 
coils used in NMR technique [ij. Such a difference 5g 
can be realized by g-factor engineering in various quan¬ 
tum dot or quantum well systems, by changing the al¬ 
loy composition in A 3 B 5 [T^ [Tg, or due to localization 
of electrons in different energetic valleys in a Ge/Si QD 
system m. A Ge/Si system with quantum dots seems 
to be more suitable, because such a pair of electrons can 
be localized at one quantum dot, one at the apex of the 
quantum dot, and a second electron at the QD base edge 

The essential difference of our approach from previ¬ 
ous schemes of quantum computation is the use of the 
g-factor difference of two electrons for effective control 


(switching on/off) of the exchange interaction between 
them. In existing works the g-factor difference is used 
for addressing individual qubits and implementation of 
single-qubit operations. There have been proposed very 
original schemes of single-qubit gate implementation, us¬ 
ing the g-factor difference, but without conventional mi¬ 
crowave pulses [3 For example, in the work of 
Levy [3 periodic modulation of the exchange coupling 
within logical qubit at the Rabi frequency D = A.gB is 
proposed to produce tt and 7 r /2 pulses. In some works 
the g-factor difference is used in spin-charge conversion 
schemes. For example, in the work of Yokoshi et al. [3 
sequential spin-to-charge conversions allow to implement 
the full Bell state measurement of electron-spin qubits. 
However, in all these proposals the exchange coupling 
must be controlled by electrodes. Our approach allows 
the implementation of single-qubit and two-qubit gates 
without electrodes controlling the exchange within one 
operational unit (an electron pair localized at one QD). 

For large-scale quantum computation at least 1000 
qubits are needed. It is possible to create a large or¬ 
dered array of quantum dots by epitaxial growth on pre¬ 
patterned substrates [2l|- Let each quantum dot possess 
a pair of electrons with a constant coupling and different 
g-factors [3 , then it can be considered as a basic ele¬ 
ment for logical gates implementation. Between quantum 
dots one can organize the SWAP operations by means 
of additional electrodes created on the surface. Since all 
quantum dots are grown in one and the same growth 
manner, they are almost identical. To provide the se¬ 
lective access to an individual QD we suggest creating 
a parallel layer of smaller QDs (a storage layer), where 
two small quantum dots correspond to one large QD in 
the first layer (Fig. 1). Let electrons localized on smaller 
QDs have the g-factor quite different from g-factors of 
electrons in the first layer, then they are always out of 
resonance, and these quantum dots can serve as storage 
elements for electron spins. In the idle mode electrons 
are located in the storage layer. The similar structures 
were created using the strain-induced nucleation of QD 
molecules above buried nanomounds [3| . For performing 
logic gates the electrons from a pair of QDs in the stor¬ 
age layer are pushed to the corresponding QD in the up¬ 
per (operation) layer by the “Address” electrode placed 
above the selected QD (see Fig. 1). The similar idea is 
described in Ref. 0 where electrons are driven by elec¬ 
trodes away from the dopant ion into layers of different 
alloy composition (and respectively different g-factors) 
for implementation of one- and two-qubit operations. 

Our approach has the following advantages: it allows 
(i) to remove a half of electrodes controlling the exchange 
interaction, namely the electrodes controlling the ex¬ 
change within electron pairs, (ii) to refuse the microwave 
pulse technics, (iii) to refuse the local strong magnetic 
field, because both magnetic fields B and SB{t) used in 
our approach are global, (iv) to use standard ESR and 
NMR techniques. 

In this work we verify the possibility of implementation 
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FIG. 1. Arrangement of quantum dots in the active (upper) 
layer and the storage (lower) layer, designed to perform quan¬ 
tum logic operations based on the control of the exchange 
interaction by means of the single spin rotations. 


Electrons are subjected to a microwave radiation circu¬ 
larly polarized in the xy-plane. The magnetic field Bu,(t) 
of the microwave is 

cosnt, sinnt, 0), 

where Bm is the amplitude of the microwave field, and If 
is the frequency. The latter is chosen to be equal to the 
mean value of Larmor frequencies of two electrons: 

nn = 


of basic logical operations on the example of the Ge/Si 
system with quantum dots. Ge QDs grown at certain 
growth conditions allow to localize two electrons on the 
same QD: near the QD apex and near the QD base edge 
(Fig. 1), that provides different g-factors due to local¬ 
ization in different A valleys, namely = 1.9995 and 
g±_ = 1.9984 [l7|. Being close to each other, these elec¬ 
trons can have a sufficient overlap of wave functions re¬ 
quired for two-qubit operations. One more advantage of 
the Ge/Si heterosystem is the location of electrons in Si, 
that can lead to long spin relaxation times due to small 
spin-orbit interaction in this material and small concen¬ 
tration of the ^®Si isotope with a nonzero nuclear spin 
[ 2 ^ [ 2 ^ . Direct measurements of electron spin lifetimes 
in dense arrays of Ge/Si QDs give the times of the or¬ 
der of 10 /iS [^. The spin lifetimes can be increased by 
switching off the main mechanism of spin relaxation in 
Ge/Si QD system—stochastic spin precession at the tun¬ 
neling between QDs [ 2 ^. Indeed, partial suppression of 
the tunneling provides fourfold increase of spin lifetimes 
[ 23 . Further increase of spin lifetimes can be reached 
by using the isotopically pure material ^®Si [ 2 ^. In gen¬ 
eral, other heterosystems, such as A 3 B 5 , can be used for 
realization of our proposal. 


II. MODEL 

Let us represent the proposed model in a mathematical 
form. We consider two electrons with a constant tunnel¬ 
ing coupling in the magnetic field (see Fig. 2). Let these 
electrons have close 5 -factors, differing on the small mag¬ 
nitude, gi= go- 5gl2 and 52 = 5o + ^5/2 with Sg go- 
A constant magnetic field Bq is applied in the z-direction, 
providing Larmor precession of the electron spins. This 
precession is described by the term 

Ho = gsBoigiSiz + g2S2z), 


thus both electrons in the stationary conditions are out 
of the spin resonance. The corresponding term in the 
Hamiltonian is 


Hjnit) = HBBjn{t){giSi + 5282 )- 


These electron spins are considered as qubits. The 
model includes the interaction between qubits having a 
form of the Heisenberg exchange interaction: 

= JS1S2, 

where J is the exchange integral, which is supposed to 
be constant. 

The main idea is to make the system controllable with 
an addition of a small time-dependent magnetic field 
SB{t). This magnetic field is also directed along the axis 
z, and the corresponding term in the Hamiltonian is 

Hc{t) = g.BSB{t){giSi^ +52-52z)- 

By means of this small additive it is possible to enter one 
of electrons into the resonance for a required time, and 
to carry out the spin turning on a desired angle through 
Rabi oscillations. Changing the additive SB{t) one can 
rotate the spin of another electron. In this way it is 
possible to manipulate with qubits in this system and to 
carry out the main quantum logical operations. 

The spin behavior becomes simpler in the reference 
frame, rotating with the frequency D. In this frame, the 
main contribution to the spin precession is removed and 
only relatively slow dynamics remains. The full Hamil¬ 
tonian of the system of two coupled electrons, written in 
the rotating frame, takes the following form: 


H{t) — Ho + Hm{0) + Hint + Hc{t) — + S 2 z)- 


where /r_B is the Bohr magneton, 51 and 52 are 5 -factors In the basis |tt)j It/-): lit): IJ-i) this Hamiltonian has 
of the electrons. Si and S 2 are their spin operators. the following matrix representation: 
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with an arbitrary phase a. Then the transformations 
themselves are related in a similar way: 

Ut=T = Ue^‘"- 

The transformation Ut=T depends on the function 
5Bit) and on the duration T of the logical operation. 
To find the best implementation of the given quantum 
gate U, we minimize the Frobenius norm of the deviation 
\\Ut=T — Ue*“|| by variation of the function 5B{t) and 
parameters T and a. 

The most suitable form of the small additive 6 B{t) is 


FIG. 2. The proposed model of the two-qubit system allowing 
the implementation of basic one- and two-qubit operations 
using a time-dependent global magnetic field B{t). 

III. REALIZATION 

In this Section we describe the strategy for finding an 
appropriate control function 6 B(t) for realization of de¬ 
sired one- or two-qubit operations. All operations are 
considered in the reference frame rotating with the fre¬ 
quency n around axis Z. 

Let us introduce an “ideal” unitary operator U that is 
defined by the desired logical operation and transforms 
an initial state \ tjjo) of the system of two spins into a final 
state U|V'o)- 

The actual evolution of the system state vector \tj}t) is 
described by the equation 

= ( 2 ) 

where the Hamiltonian H{t) depends on 6 B{t) according 
to Eq.(l). An unitary operator Ut, that transforms \tjjo) 
into IV't), is defined by the equation 

= H{t)Uu (3) 

and the initial condition that Uq is the identity operator. 

At some time T, the transformation Ut=T should be¬ 
come equivalent to the “ideal” operator U, in order to 
perform desired logical operation without an error. In 
other words, for all \4’) the results of actions of Ut=T and 
U should be physically equivalent: 

C/t=T|^) =U|V') 


5B{t) = Acos{ijjt + (f) + C (4) 

with parameters A, w, C depending on the type of the 
desired logical operation. Such dependence can be easily 
realized experimentally. Changing these parameters, one 
can perform various one- and two-qubit quantum gates, 
as it will be demonstrated below. 

Fig. 3 demonstrates the example of behavior of error 
functional 

/(t) = min„||t/t-Ue*“f, (5) 

during the operation “•7r/2-rotation around Z". It is 
clearly seen that the error functional has the minimum 
value at t = T. 

One can also use a more complicated modulated func¬ 
tion SB{t) having zero values at the beginning and at the 
end of the control action: 

dBit) = Asin(woO(cos(a;t -\- tp) + C), (6) 

where wq = 2tt/T with T being the time of gate imple¬ 
mentation. Such a form of the control function is con¬ 
venient from the standpoint of the implementation of a 
large sequence of logical operations and smooth interfac¬ 
ing between operations. The search for optimal parame¬ 
ters (wq, a, w, p, C) of the function (6) can be done in 
two stages. The first step is the solution of the problem 
with a control function described by Eq. (4) to determine 
the time T at which the error of gate implementation is 
minimal. At the second step one can apply the func¬ 
tion described by Eq. (6) with wq = 27r/T to refine the 
solution. 

In the following Section we will demonstrate that at 
reasonable values of tunneling coupling and g-factor dif¬ 
ference one can realize a minimal set of one- and two- 
qubit operations required for quantum computation. As 
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FIG. 3. Time dependence of the error functional f{t), Eq.(5), 
for the operation “7r/2-rotation around Z”. The time, when 
the error functional has the minimum, is indicated as T. Inset 
shows the tolerance range for the time T. One can see that the 
end of logical gate can be controlled with accuracy ±2 ns, at 
that the increase of error / does not exceed threshold Sfmin- 
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we mentioned above, the small difference of g-factors 
dg ^ 10“^ observed recently in Ge/Si system with quan¬ 
tum dots [a is suitable for our purposes. The mag¬ 
nitude of the exchange interaction can be estimated as 
J ^ AP/U ~ 10“^° eV, where / is the tunneling inte¬ 
gral, I ^ 10“^ meV, and U is the Coulomb interaction, 
U ^ 30 meV. The microwave magnetic field can be taken 
as Bjn = 1 G. The time-dependent magnetic fields create 
the electric fields, but their magnitude is rather small. 
The alternating magnetic field with an amplitude ^ 1 G 
causes the electric field of the order of 300 V/cm. On the 
quantum dot length scale 10 nm) the induced voltage 
is about 0.3 mV. Compared with the electron confine¬ 
ment potential in QD system 100 meV), this value is 
negligible, and we do not take its effect into account. 

IV. RESULTS AND DISCUSSION 

In Table I the time T and the value of gate fidelity 
F are given for all verified logical operations. The gate 
fidelity is the probability of obtaining the correct result, 
averaged over all initial states [^. The connection 
between the error functional / and fidelity F is described 
in Appendix A. For all realizations we use the function 
dB{t) described by Eq.(4). Parameters A, uj, tp, C were 
chosen to minimize the deviation of the fidelity from 1, 
using the optimization methods described in Ref. [12. 

Fig. 4 demonstrates the results for the operation 
“quantum information storage”. The top panel shows 
the evolution of spin components Sx, Sy, Sz for the 
first electron. For the second electron the spin evolution 
is shown at the central panel. And the bottom panel 


FIG. 4. Time evolution of spin components of the first (top 
panel) and the second electron (center panel) during “Stor¬ 
age operation”. The blue solid curves are related to the 
^-components of the electron spins Si and S2, the green 
dashed curves are related to the S'^-components, the red dot¬ 
ted curves are related to the Sz-components. The control 
magnetic field 5B{t) is shown in the bottom panel. 


demonstrates the time dependence of the magnetic field 
SB{t), that provides the conservation of the system in 
the initial state. It is clearly seen that the system un¬ 
der the action of the controlling magnetic field returns 
to its initial state. The different behavior is observed for 
the first and the second electrons, the first electron be¬ 
sides high-frequency oscillations of spin components has 
low-frequency ones, while the second electron has only 
high-frequency oscillations. This difference is explained 
by nonzero constant C in the function 6 B{t), that shifts 
one of the electrons closer to the resonance. In principle, 
it is possible to realize the identical behavior of the spins 
with C = 0, certainly in this case the time T and other 
parameters of the function SB{t) will be changed. 

In Fig. 5 the evolution of the electron spins without 
time-dependent controlling magnetic field is shown. Re¬ 
sults demonstrate that it is impossible to find such a mo¬ 
ment of time, at which both electron spins return to their 
initial orientations. 

The evolution of electron spins for the operation “'k/2- 
rotation around Z” is presented in Fig. 6 for the follow¬ 
ing initial conditions. The first electron spin is oriented 
along the axis X, the second electron spin is directed 
along axis Z. As a result of the operation, the spin of 
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TABLE I. Results of numerical experiments. The controlling magnetic field has the form 5B{t) = Acos{u}t -\- ip) -\- C. The 
value (1 — E) gives the averaged probability of the error, F is the gate fidelity. The results were obtained for the following 
parameters: mean g-factor value go = 2, g-factor difference 5g = 1.1 ■ 10“®, exchange interaction J = 10“^° eV, microwave 
frequency fl = 9 ■ 10® Hz, microwave field amplitude Bm = 1 G. 


Operation 

1 - 

F 

u (MHz) 

A(G) 

V 

(■ Ml, 

r(/rs) 

Storage 

1.95 • 

lo-'^ 

14.80 

1.24 

0.08 

2.417 

7.212 

VSWAP 

9.44 • 

10"'‘ 

6.79 

0.99 

3.14 

-0.009 

11.104 

SWAP 

1.06 • 


6.79 

0.87 

0.00 

-0.005 

22.209 

7r/2-rotation around Z 

4.00 • 

10"® 

14.39 

0.81 

0.85 

2.319 

1.700 

7r/2-rotation around X 

3.64- 

10”® 

17.43 

1.96 

6.27 

-2.710 

9.375 

7r/2-rotation around Y 

4.53 • 


14.40 

1.24 

1.82 

1.323 

1.309 

7r/4-rotation around X 

1.40 • 

10'® 

9.741 

0.07 

3.05 

0.253 

0.664 

7r/8-rotation around X 

5.55 • 

10”® 

4.816 

1.00 

5.99 

-0.029 

1.424 
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FIG. 5. Free evolution of spin components (SB = 0) of the 
first (top panel) and the second electron (bottom panel). The 
system does not return to the initial state. 


the first electron has arrived to the orientation along the 
axis V, while the second electron spin has returned to its 
initial orientation. This demonstrates the correctness of 
the proposed gate implementation. 

Fig. 7 demonstrates the system evolution for the 
SWAP operation. This operation is implemented if the 
initial state transforms to the final state in the following 
manner: 


Itt) ^ Itt) e“, 

in) ^ in) e“, 
lit) ^ Iti) e“, 

lit) —>• lit) e*“, 

where horizontal arrows represent the transformation 
from m) to nr) under the given Hamiltonian. In Fig. 7 
the hrst electron was initially in the spin-down state, and 
the second one was in the spin-up state. As result of oper¬ 
ation the electrons have exchanged their spin directions. 




time (ps) 

FIG. 6. Time evolution of the spin components of the first 
(top panel) and the second electron (center panel) during the 
operation “7r/2-rotation around Z". The blue solid curves 
are related to S'a;-components of electron spins Si and S 2 , 
the green dashed curves are related to S'^-components, the 
red dotted curves are related to 5 ' 2 -components. The control 
magnetic field 5B(t) is shown in the bottom panel. 


Simple reduction of the operatio n time b y a factor of 
2, T = Tswap/‘2', produces the VSWAP gate, which 
has then the maximum entangling capability. The re¬ 
sults presented in Table I show that the frequency w is 
the same for the V SWAP and SWAP gates, when the 
amplitude A and the phase tp are slightly different. 

The error probability (l — F) for the spin rotation gates 
is found in the range (1.4 -f 5.5)-10“^, while the error 
probability for the SWAP operation is approximately 






































































7 





time(ns) 


FIG. 7. Time evolution of Sz-components of the first (top 
panel) and the second electron (center panel) during the 
SWAP operation. The control magnetic field 5B{t) is shown 
in the bottom panel. 


10“^. We suppose that the main source of errors in single 
qubit operations (spin rotations) is a closeness of electron 
^-factors. To verify this hypothesis we perform numerical 
experiments for the operation “ 7 r/ 2 -rotation around X” 
with twofold increased <5^ = 2.2 -10“^. As a result we 
obtain the five-fold reduction of the error, (1 —F) ~ 10“^. 

A quality of the SWAP operation depends on two pa¬ 
rameters: exchange interaction J and detuning 5g/2 of 
the electron g-factors from the resonance. If the magni¬ 
tude of the exchange will be very large, the evolution of 
the two-electron system becomes less controllable. For 
example, the numerical experiment for the SWAP gate 
with J = 10“® eV (at 5g = 1.1 • 10“^) gives the error 
{I — F) = 3.52 • 10“^. Regarding the parameter 6 g, its 
increase also leads to the error rise. In particular, the 
SWAP operation with 6 g = 2.2 • 10“^ (J = 10“^° eV) 
has the error {1 — F) = 1.73 • 10“^. In the last case the 
detuning of electrons from resonance is too strong and 
it becomes hard to use the spin rotations for control (or 
elimination) of the exchange interaction. In Fig. 9 the 
behavior of the error (1 — F) with change of Sg and J is 
shown for the SWAP operation. It is clearly visible that 
the values of 6 g = 1.1 • 10“^ and J = 10“^° eV used in 
our numerical experiments are close to the optimal ones. 
It should be noted that we chose these values based on 
the parameters of real Ge/Si quantum dots. From the 
dependence shown in the left panel of Fig. 9 it seems 



C(G) 


FIG. 8 . Dependence of the error functional / on the deviation 
of the parameter C (a constant shift in the controlling field 
SB(t) = Acos{ujt + tp) C) for the '/SWAP gate. Inset 
shows that at AC < 0.01 G the error / does not exceed the 
threshold Sfmin- 
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FIG. 9. Dependence of the error (1 — F) on the g-factor 
difference 5g (left panel) and on the exchange integral J (right 
panel) for the SWAP gate. Here F is the gate fidelity. 


that the value Sg = 5 ■ 10“'* is more appropriate than 
Sg = 1.1 • 10“^. However, the decrease of Sg complicates 
the addressing individual qubits, that provokes the rise 
of the error (1 — F) during single-qubit operations. Thus, 
the search for the optimal values of Sg and J should be 
carried out with respect to the whole set of one- and 
two-qubit operations. 

It is worth noting that each optimization problem con¬ 
sidered in this paper can have a few possible solutions 
with approximately equal errors. In principle one can 
use any of these solutions. Each solution has its own 
unique moment of time, when the system reaches the de¬ 
sired state. Since the behavior of the error functional / 
represents very strong oscillations, the success of our ap¬ 
proach depends on the accuracy of time determination. 
The tolerable error At, at that the averaged probabil¬ 
ity of error (1 — F) stays the same order of magnitude, is 
about ±2 ns (Fig. 3). Therefore the gate duration should 
be controlled with accuracy ~ 1 ns. Such a requirement 
can be easily satisfied with modern experimental equip- 
















































































































merit. 

Additionally we have studied the reliability of our ap¬ 
proach with respect to deviation of parameters of the con¬ 
trolling field 5B{t) {A, w, (/>, C) from their optimal values. 
The study was performed for the V SWAP operation. In 
a real experiment all these parameters can be controlled 
with some accuracy, and we have found for each parame¬ 
ter the tolerance range at which the fidelity F remains of 
the same order as the optimal value. We define the toler¬ 
ance range from the condition / < where fmin is 

the error functional value given in Table I. In Fig. 8 the 
tolerance range for the parameter C is shown. From the 
inset to this figure one can see that the tolerance range 
AC = ±0.01 G. However, the fidelity is more sensitive to 
the amplitude A. The deviation AA = ±0.001 G causes 
the increase of the error / up to 3fmin- The frequency 
u! must be controlled with accuracy Aw = ±5 kHz. The 
tolerant change of the phase (p is A(p = ±0.12. It should 
be noted that within the indicated tolerance ranges the 
time of the gate implementation does not change. 

Finally we would like to discuss the experimental feasi¬ 
bility of the magnetic field control with required accuracy. 
The magnetic field in our proposal consists of two parts 
B ± 6 B{t). The first large part is of the order of 10^ G, 
the second small part is about 1 G. Both fields should 
be stabilized with accuracy 10“^ G. The most challeng¬ 
ing problem is stabilization of the large magnetic field, 
because the relative accuracy should be of the order of 
10“®. Modern superconducting magnets are able to meet 
this requirement. For example, the specification of uni¬ 
formity of the magnetic fields in the general NMR appa¬ 
ratus is O.OI ppm in a sample space and 0.01 ppm/hour 
in time-based stability. The additive part 6 B{t) can be 
produced by a small extra coil, and should be controlled 
with a reasonable accuracy ~ 0 . 1 %. 

V. CONCLUSIONS 

We have proposed a new method of the implementa¬ 
tion of one- and two-qubit gates in a system of two elec¬ 
trons with constant exchange coupling. The only con¬ 
trolled parameter is the small time-dependent addition 
to the magnetic field. Small g-factor difference between 
two electrons provides not only the selective access to the 
individual qubits, but also the effective control of the ex¬ 
change interaction between them. We have verified the 
possibility of basic one- and two-qubit operations in the 
frame of proposed model and found the parameters of 
controlling magnetic field allowing to perform these op¬ 
erations with an acceptable accuracy. 
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Appendix A: relation between error functional and 
fidelity 

In this Appendix, we will show how to find the gate 
fidelity F from the functional / defined by Eq. (5). 
Frobenius norm of a square matrix A is defined as 

||A|| = v/Tr(AA+), 

therefore 

||C/t - Ue“f = Tr [(Gt - Ue“) (t/+ - U+ e"”)] 

= Tr {UtU+) + Tr (UU+) - 2 Re [e“ Tr (C/+U)] . 

In the latter expression, UtU^ and UU’*' are identity ma¬ 
trices 4 X 4, so their traces are equal to 4. Denoting the 
remaining trace as Z, 

Z = Tr (t/+U) , 

one obtains the following relation: 

||C/t - Ue*“f = 8 - 2Re (e“ Z) . 

Hence, according to Eq. (5), 

f{t) = min„ \\Ut - Ue*“f =8-2 jZj. (Al) 

Since the matrix U^U is unitary, its four eigenvalues 
can be expressed as 6*“^% ..., where pi,. ■ ■ ,^4 are 
real numbers. So the trace Z of this matrix is equal to 

4 

Z = ^F^\ 

therefore 

k I 

= 4 ± 2^cos((^fc - </J/) = 4 ± 2M, (A2) 

k<l 

where 

M = ^cos{(pk - (fi)- 

k<l 

Using Eqs. (Al) and (A2), one can express the quantity 
/ via M: 

/ = 8 - 2V4 ± 2M. (A3) 

Now let us express the gate fidelity F via M. The gate 
fidelity is the probability, averaged over all initial state 
vectors \ip), that the system will be found in the desired 
state U|^/’) after passing through the gate: 


F = \{i;\U+m;)\ , 
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where the overline denotes averaging over initial states. 
Let us consider this expression in a basis that diagonalizes 
the matrix 


C/+U = 


therefore 

= X! + 2^ |cfe|2|ciPcos((/5fc - Lpi). 


k I 


k k<l 

It is possible to show that 


/ 

0 

0 

0 

\ 


/ 

Cl 
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[l/lO 

if k 
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0 

0 




C2 



0 

1 — 1 

II 

u 

1 

if k 

0 

0 


0 


- IV’) = 


C3 


5 












therefore 



0 

0 

0 

g*¥>4 

J 


1 

C4 

J 





F = 


4 + M 
10 


(A4) 


Comparing Eqs. (A3) and (A4), one can easily see that 


where ci,..., C 4 are complex amplitudes that define the 

state !'(/')■ In this basis, If f ig small enough, then (1 — F) is approximately pro¬ 

portional to /: 


(V^|C/+U|V^)=^|cfcpe*^^ 

k=l 


Hence, the minimization of the error functional / is 
equivalent to the minimization of the deviation of the 
fidelity F from 1. 
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